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We investigated phase defects in a quasi-one-dimensional commensurate charge density wave 
(CDW) system, an In atomic wire array on Si(lll), using low temperature scanning tunneling 
microscopy. The unique four-fold degeneracy of the CDW state leads to various phase defects, among 
which intrinsic solitons are clearly distinguished. The solitons exhibit a characteristic variation of the 
CDW amplitude with a coherence length of about 4 nm, as expected from the electronic structure, 
and a localized electronic state within the CDW gap. While most of the observed solitons are 
trapped by extrinsic defects, moving solitons are also identified and their novel interaction with 
extrinsic defects is disclosed. 



A topological soliton is a local solitary wave bridging 
two degenerate states, which has been an important con- 
ceptual tool in many branches of science pQ . For exam- 
ple, in electronic systems, solitons are responsible for the 
high electric conductivity in conjugated polymers such 
as polyacetylene [HE]. A soliton exists as a nontriv- 
ial phase defect separating energetically degenerate one- 
dimensional (ID) charge density wave (CDW) states with 
an electronic state inside the CDW gap [JHS]- While 
solitons in CDW systems show up their existence mainly 
through their consequences in transport and optical prop- 
erties, it has been challenging to observe them directly 
in real space mainly because of their microscopic dimen- 
sion and high mobility. This obstacle has prevented ex- 
perimental investigations on various microscopic inter- 
actions of solitons predicted theoretically, such as those 
with other solitons and defects [7Vll2j. 

The experimental difficulty can partly be overcome 
by the high spatial resolution of scanning tunneling mi- 
croscopy (STM). Only very recently, a soliton was visu- 
alized by STM on an incommensurate ID CDW system 
of NbSe3 [13j . which was immobilized by an unknown 
reason. On the other hand, for a commensurate CDW 
system, a soliton has not been clearly identified partly be- 
cause it can easily be confused with non-solitonic phase 
defects such as structural defects. Indeed, recent STM 
studies of metallic atomic wires on silicon surfaces in their 
CDW ground states revealed various phase defects [TH - 
[T5] , However, the origin and solitonic characteristics of 
these phase defects have been debated and largely un- 
clear [HI ED]. 

In this Letter, we carefully reinvestigated with STM 
various local phase defects of In atomic wires on Si(lll) 
in their quasi-lD CDW states PHM ITS] . We examine 
the solitonic characteristics of these defects such as their 
shapes (local variations of CDW envelopes), electronic 
states, and mobility. We can unambiguously distinguish 
topological solitons from various non-solitonic phase de- 
fects, solving the current debate on this system [191120]. 
This opens up the possibility of investigating microscopic 



interactions of solitons in a commensurate CDW system 
and, in particular, the present work elucidates novel in- 
teractions of solitons with extrinsic phase defects. 

The experiment was carried out with a ultrahigh- 
vacuum cryogenic STM (Unisoku, Japan). The 
Si(lll)4x 1-In surface of In atomic wires was prepared 
by depositing 1 ML of In onto the clean Si(lll)7x7 sur- 
face at an elevated temperature |21j . Subsequently, the 
sample was cooled down to 78 K, well below the CDW 
transition temperature of ~ 125 K, for STM measure- 
ments. All STM images presented here were taken in the 
constant-current mode with typically a tunneling current 
of 50 pA and a sample bias voltage of —1 V. 

An In wire of Si(lll)4x 1-In has two In atomic chains 
and is separated by Si chains with periods of 4ao and lag 
perpendicular to and along the wire, respectively (do is 
the Si lattice spacing, 0.384 nm) [21] [22] . Upon cooling, 
a commensurate 4x2 CDW state develops through the 
periodicity-doubling distortion along a wire and shows a 
clear gap opening for a half-filled metallic band [33]. In 
addition, interwire coupling forces CDW wires to have 
opposite CDW orientations alternatively perpendicular 
to the wire, leading to 8x2 ordering [55]. Since there 
arc four different ways to form a 4x2 CDW state from 
two In atomic chains within a single wire [see the inset of 
Fig. 1(a)] j2~4] . we can expect up to three different kinds 
of phase defects (and their mirror symmetric versions) as 
described in Fig. 1(a). One is a phase slip defect (PSD) 
between two-fold degenerate CDW states with the same 
CDW orientation, which are shifted by ir or half a CDW 
period (loo) along the wire [left panel in Fig. 1(b)]. On 
the other hand, two CDW states with opposite CDW 
orientations meet at a phase flip defect (PFD). This de- 
fect is unconventional but comes from the unique double 
chain structure of the In wire; only one of two In atomic 
chains contains a ir phase shift. Unlike the PSD, the 
PFD lacks the two-fold degeneracy under the interwire 
coupling. Thus a PFD would induce two long domain 
walls [dashed lines of middle panel in Fig. 1(b)] along the 
wire. Instead, the PFDs can exist near the edges of CDW 
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FIG. 1. (a) Schematic diagrams of three possible phase de- 
fects (N = 1) on a 4x2 CDW wire (solid ovals). Inset shows 
four degenerate CDW states in different colors for a single 
wire with two In atomic chains (black dots). The atomic 
structure is largely simplified [23] . (b) Schematics of possible 
phase defects: a phase slip defect (PSD, left), a phase flip 
defect (PFD) with the broken 8x2 order (middle), and a line 
boundary of PFDs (right). Two red ovals indicate two out of 
four possible 8x2 configurations at the given CDW orienta- 
tion. STM images showing (c) two PSDs and (d) PFDs and 
composite defects with (e) a corresponding differential con- 
ductance map taken at +0.15 V |22j . Two nearest normal 
CDW maxima are used to determine lengths of defects as de- 
noted by slashes. The distinct common local features of the 
defects are marked by dots [(c)-(e)]. 



wires such as step edges or cluster into a line boundary 
[right panel in Fig. 1(b)] of two CDW domains across the 
wires while various isolated PSDs are observed [22j . 

Figures 1(c) and 1(d) show STM images of various 
phase defects, which are static enough to be imaged by 
STM. Two PSDs in Fig. 1(c) show contrasting length 
scales. A short PSD has a distinct local structure as in- 
dicated by dots [Fig. 1(c)], which abruptly induces the 
CDW phase shift. In contrast, a long PSD in the neigh- 




FIG. 2. STM line profiles (dots) of (a) a long PSD and (c) a 
long PFD obtained along the dashed lines in Figs. 1(c) and 
1(d) together with their best fits by the soliton shape with 
(red) and without (green) considering the lattice modulation. 
The line profiles are high-pass filtered to get rid of slowly vary- 
ing backgrounds. With the lattice modulation, the soliton 
shape is modified as Zg(x) = Z${x) ± Al cos[27r(a; — s5o)/ao], 
where Al is the height modulation of underlying atoms. The 
inclusion of the lattice modulation results in a much better 
fit with almost the same coherence length. £( a ) = 3.75 ± 0.14 
(3.75 ±0.17) nm and £ (c) = 3.50 ±0.23 (3.40 ±0.26) nm with 
(without) the lattice modulation, (b) and (d) Enlargements 
of (a) and (b) around the soliton center. 



boring wire gradually changes the CDW phase over a 
relatively long distance of ~ 20a without such a dis- 
tinct local structure. The latter length is consistent for 
most of the long PSDs observed. On the other hand, as 
mentioned above, the PFDs tend to form a line bound- 
ary perpendicular to the wires [Fig. 1(d)] [22]. The short 
PFD has a similar local structure [dots in Fig. 1(d)] to 
a short PSD. A long PFD changes its phase gradually 
like a long PSD but for only one In chain in the wire. 
The short PFDs as well as the short PSDs seem to act as 
building blocks to form various composite phase defects 
as shown in the top wire of Fig. 1(d). 

The long PSDs are similar to what was assigned as a 
soliton by Morikawa et al. [2] and the short PFDs were 
argued as solitons by Zhang et al. mainly due to the phase 
shift involved [T8H2Pj] . Note, however, that various ex- 
trinsic (structural, chemical, and charge) defects can in 
principle cause CDW phase shifts, which should be dis- 
tinguished from a intrinsic soliton as the lowest-energy 
excitation from the CDW ground state [13] ■ In terms of 
the shape, such a soliton excitation would have a well de- 
fined shape Zs(x) — Z a ± Ac tanh[(x — Xo)/£] sin[7r(x — 
Xo)/ao\, where Ac is the CDW amplitude and £ is the 
microscopic coherence length [4] [13] . The CDW ampli- 
tude varies smoothly within the coherence length, half 
the apparent size of a soliton, which is determined by 
£ = 7ui.f/A, where Vf is the Fermi velocity and A is half 
the CDW gap j3] [25] ■ For the present case, the band 
dispersions [23] and the gap size [22] yield £ = 3.4 nm 
~ 9ao, which is similar to that of the soliton in poly- 
acetylene, ru 7a It is thus apparent that the short 
PSDs and short PFDs are not solitons judging from their 
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shapes and lengths (< 2£~18ao). They are most likely 
structural defects or adsorbate-induced ones. Structural 
defects with a ir phase shift in the dimerized chain struc- 
ture were observed in various systems including the clean 
Si(001)2xl [35] and Si(lll)5x2-Au surfaces [27]. 

However, the long phase defects are clearly differ- 
ent. We fitted a line profile obtained from a long PSD 
[Fig. 1(c)] with the above soliton shape. As shown in 
Fig. 2(a), the profile is excellently fitted; the gradual de- 
cay of the CDW amplitude and the ir phase shift around 
the center are clearly visualized. The coherence length 
fitted is £ = 3.75 ±0.17 nm [Fig. 2(a)] in good agreement 
with that estimated above. This form is very much con- 
sistent with the trapped soliton observed in NbSe3 [13]. 
We thus interpret that a long PSD is an intrinsic soliton 
immobilized by some trapping potential. The trapping 
potential is obviously provided by a neighboring short 
PSD [indicated by the wavy line in Fig. 1(c)] since most 
of long PSDs are observed as paired with neighboring de- 
fects. This trapping potential affects the soliton shape to 
deviate marginally from its theoretical form as indicated 
by the ovals in Fig. 2(a). The deformation is more sig- 
nificant on the In atomic chain closer to the neighboring 
defect [22]. 

A profile of a long PFD can also be consistently fitted 
with the same function with a similar coherence length 
[Fig. 2(c)] [22] ■ It suggests that the long PFDs are also 
similar solitons. However, the deviation near a neighbor- 
ing defect, which is located near the center, is apparently 
larger in this case [Fig. 2(d)]. This may be explained by 
a CDW perturbation due possibly to the Friedel oscilla- 
tions induced by a neighboring defect [the wavy line in 
Fig. 1(d)] as also observed in NbSe3 [13] . Including the 
local electronic response to the neighboring defect, we 
could reproduce precisely the shape of a PFD with the 
soliton profile [22] . 

Since a soliton would have an electronic state within 
the CDW gap [TJHHS], we checked this localized electronic 
state by scanning tunneling spectroscopy for short and 
long PFDs [22 . From the differential tunneling conduc- 
tance map obtained just below the upper CDW gap edge 
[Fig. 1(e)], we clearly found that the long PFD shows a 
localized state (~ +0.15 eV) within the band gap with 
a similar length to the coherence length in contrast to 
the short PFDs and their composite defects 22 . The 
deviation of this soliton state from the Fermi energy may 
be ascribed to the structural or electronic perturbation 
on the long PFD by the neighboring defect or a strong 
interaction with the underlying lattice in a commensu- 
rate CDW system, where the CDW phase is completely 
locked to the underlying lattice. On the other hand, a 
similar spectroscopy measurement on a long PSD was 
hardly possible since it can rather easily be detrapped 
during a long measurement time for spectroscopy. 

The trapping or detrapping assumes another impor- 
tant characteristic of a soliton to be addressed, the mobil- 




FIG. 3. (a)-(c) A series of STM images taken consecutively 
with their start times indicated. Two defects indicated by 
white and black triangles are a reference and a hopping PSD, 
respectively. Sometimes, two abrupt CDW phase shifts (0 — > 
7r — > 2n) occur in a short time span as marked by a black 
arrow in the inset of (c). (d) Stacked STM line profiles (the 
newest line profile at top most) taken along a dashed line 
in (c). Each black arrow indicates the two successive CDW 
phase shifts. 



ity. Therefore, the assignment of a trapped soliton would 
not be complete without observing freely-moving soli- 
tons. In order to find any moving solitons, we scanned the 
same restricted area repeatedly. Figures 3(a)-(c) are a 
series of sequential STM images showing the same CDW 
wires around two composite defects with (a 15ao-PSD) 
and without (12ao) a phase shift. One can observe that 
the 15ao-PSD hops by 2a to the right [Fig. 3(b)] and 
then hops back to the original position in a time scale 
of a few minutes [Fig. 3(c)]. Most of the short defects 
showed a similar random hopping motion by a multiple 
of a CDW period (2ao) at a time. This motion is quite 
similar to that previously reported by Zhang and oth- 
ers [18] . However, such a hopping motion is natural for 
trivial defects like the non-solitonic short PSDs. 

In sharp contrast to the low frequency hopping, we 
could observe a much faster (beyond our temporal res- 
olution) shift of a CDW phase as marked by the black 
arrow in the inset of Fig. 3(c). A ir phase shift occurs 
abruptly but the phase returns back in a short time of 
~ 50 msec. This abrupt phase shift is more clearly visu- 
alized by continuously scanning the same line along one 
CDW wire. Figure 3(d) shows an example of such a line 
scan series, where two events of abrupt phase shifts are 
observed (indicated by black arrows). Each event consists 
of two successive ir phase shifts (0 — > ir — > 2ir). These 
fast phase shifts may result from a moving soliton that 
moves too fast to be captured with STM. Then, the above 
STM image would indicate two solitons passing sequen- 
tially the scanning tip. Rather frequently, these abrupt 
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FIG. 4. (a) Stacked STM line profiles showing the hopping 
of the PSD. Dotted lines and white arrows indicate the posi- 
tions of the PSD and the hopping directions, respectively, (b) 
The STM profiles extracted from (a) are displayed in a time 
sequence, with the earliest on the top (0) and the latest on 
the bottom (11). Left (right) curves show the line profiles ob- 
tained by scanning from left (right) to right (left). The scan 
speed was set to 150 msec per line. Black triangles and verti- 
cal arrows indicate the positions of the PSD and abrupt CDW 
phase shifts, respectively. Gray boxes highlight emerging long 
PSDs or solitons. (c) Fitting of a soliton's profile (dots) with a 
soliton solution (solid line, £ = 3.56±0.48 nm). Note that raw 
data was used for fitting in contrast to the trapped solitons. 
(d) Schematics of the interaction between solitons (51,2) and 
a PSD. Black and white areas indicate different CDW states. 



see the scan (1) and (2)] and disappeared after another 
7r phase shift [ir —t 2ir, (3) and (4)]. This indicates that 
these PSDs are moving solitons trapped transiently and 
the abrupt CDW phase shifts are created when these soli- 
tons are trapped or detrapped. 

The transient trapping of a moving soliton can be ex- 
plained by the soliton-barrier model. Theoretical simula- 
tions show that the scattering of a soliton with a poten- 
tial barrier is nearly elastic [9H11] . For an initial velocity 
Vi smaller (larger) than the critical velocity v c , a soliton 
reflects back (transmits over the barrier) [Fig. 3(d)]. In- 
terestingly, at Vi ~ v c , a soliton interacts with the barrier 
slowly and stays for a while near the barrier until it fi- 
nally escapes. If the present soliton at a velocity close to 
v c interacts with the potential barrier of a short PSD, the 
observed motion can be explained qualitatively well. One 
thing not captured in this theory is that in our case the 
potential barrier itself moves when a soliton transmits 
through it [Fig. 4(d)]. It is natural that a soliton with 
a 7T phase shift in itself translates a static phase defect 
by 7r or la . Since the soliton motion and the hopping 
of extrinsic defects are coupled, it might require an extra 
energy cost for a soliton transmission. 

It is noteworthy that most of the hopping events ob- 
served are 2ao-hoppings consisting of two lao-hoppings 
in a very short time span [Fig. 4(d)]. This suggests that 
two solitons, that is, a soliton and an antisoliton, are 
combined to move as predicted in theory [T]. This needs 
further investigation along with the issue of separating 
pure thermal hoppings of the defects from the soliton- 
induced ones. 
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CDW phase shifts occur together with the hopping of the 
PSD as shown in Fig. 4(a). For a more detailed analy- 
sis of these complex motions, we extract individual line 
profiles from Fig. 4(a), especially, at close range of the 
hopping of the PSD. We found that the PSD (black trian- 
gles) hopped by lag right after an abrupt n phase shift 
(black arrows) and two successive phase shifts induced 
the PSD to hop by 2a in total between the scan (0) and 
(4) [Fig. 4(b)]. The same successive but reverse hoppings 
occurred between the scan (7) and (11) with two abrupt 
phase shifts. 

More interestingly, we observed a few long PSDs tran- 
siently [indicated by gray boxes in Fig. 4(b)]. The de- 
tailed line profile [Fig. 4(c)] clearly indicates that they 
are the same as the static solitons discussed above. As a 
rule, they appeared after an abrupt tt phase shift [0 — > 7r, 
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